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Abstract 

'66 ' 

O ' We derive an energy density functional for non-relativistic spin one-half fermions in the limit of a diver- 

o, 

CN \ gent two-body scattering length. Using an epsilon expansion around d = A — z spatial dimensions we com- 

Q ' pute the coefficient of the leading correction beyond the local density approximation (LDA). In the case of 

^^ ■ A'^ fermionic atoms trapped in a harmonic potential this correction has the form E = Eida{^ +Cs{3N)^'^I^), 

where ^lda is the total energy in LDA approximation. At next-to-leading order in the epsilon expansion we 

find Cs = 1.68, which is significantly larger than the result for non-interacting fermions, Cs = 0.5. 
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I. INTRODUCTION 

Density functional theory (DFT) |llil2|] is widely used in condensed matter physics and quantum 
chemistry to treat quantum many-body systems. DFT relies on the fact that one can demonstrate 
the existence of an energy functional that depends only on the fermion density, and not on the 
wave functions, that has the property that the ground state energy can be obtained by minimizing 
the functional with respect to the density. There is renewed interest among nuclear theorists in 
developing a "Universal Nuclear Energy Density Functional" [3] which describes nuclear proper- 
ties all across the isotopic chart, including very neutron rich nuclei far from the valley of stability. 
Ideally, this functional should be derived from a systematic theory of the nucleon-nucleon interac- 
tion. A modern approach to nuclear forces is provided by the effective field theory (EFT) method 
ll4tl5il6|,LZ|]. EFT starts from the most general local lagrangian containing nucleon and pion fields 
that respects the basic symmetries of QCD. A power counting scheme determines the number of 
terms in the lagrangian, as well as the number of diagrams, that have to be kept at any given order 
in a low energy expansion. 

If the effective field theory is perturbative then there are systematic methods for determining 
the energy density functional [|8|, |9|]. However these methods are very cumbersome if long-range 
forces, such as pion exchange or collective modes, are important. In the case of non-perturbative 
EFT's no systematic calculations of the energy density functional exist. In nuclear matter, and in 
cold fermionic gases near a Feshbach resonance, an important source of non-perturbative physics 
is the large two-body scattering length. 

In the present work we try to address some of these problems by studying the energy density 
functional for a dilute system of non-relativistic spin 1/2 fermions with an infinite two-body scat- 
tering length. Because the s-wave cross section saturates the unitarity bound this limit is often 
referred to as the "unitarity limit". The energy density functional in this limit is important for the 
study of neutron star crusts and neutron halos in nuclei. It can also be used to describe trapped 
fermionic atoms in the vicinity of a Feshbach resonance. The Fermi gas in the unitarity limit 
exhibits a number of interesting non-perturbative phenomena. It is a superfluid, and the ratio of 
the gap over the Fermi energy is large. Superfluidity implies that the U{1) phase symmetry is 
spontaneously broken and the low energy or momentum response is carried by Goldstone modes. 

We shall compute the energy density functional up to next-to-leading order (NLO) in an expan- 
sion in derivatives of the density. Our procedure is based on an effective lagrangian for the Fermi 



m- 



gas in the unitarity limit derived in [llOn . We will determine the coefficients in this lagrangian using 



an epsilon expansion around d = 4 — e spatial dimensions [IllL [120. As a by-product we compute 
the phonon dispersion relation and the static susceptibility at NLO in the epsilon expansion. Our 
result for the energy density functional is rigorous when applied to infinite systems in which the 
density varies smoothly, but there are some limitations in the case of finite systems with a sharp 
surface, such as fermions confined in a harmonic trap. Pairing leads to an odd-even effect in the 
dependence of the energy on the number of particles, which a local energy density functional that 
depends only on the particle density cannot describe. Also, the gradient expansion breaks down 
near the surface of the system and the expansion of the energvin inverse fractional powers of the 
number of particles cannot be pushed to arbitrarily high order lIlQII . We shall discuss some possible 
approaches to overcome these limitations in Sec.lVl 

II. EFFECTIVE LAGRANGIAN AND ENERGY DENSITY FUNCTIONAL 

The energy density functional describes the response of the system to smooth variations in 
the density. This functional can be related to the effective lagrangian that governs the response 
to slowly varying external fields. The effective Lagrangian at NLO in derivatives of the external 



potential is [llOh 



£ = com3/2z5/2 + cimi/2i^ + ^ [(v2(p)2-9mVVl Vx , (1) 

VX y/in\- J 

where we have defined 

(V(p)2 

2m 
The lagrangian contains the Goldstone boson (phonon) field cp(-?,?), the chemical potentially, and 
the external potential V{x,t). The mass of the fermion is denoted by m. The functional form 
of the effective lagrangian is fixed by the symmetries of the problem, Galilean invariance, t/(l) 
symmetry, and conformal symmetry. The NLO effective lagrangian is characterized by three di- 
mensionless parameters, co,ci,C2. These parameters can be related to physical properties of the 
system. The first parameter, co, can be related to the equation of state. We have 

25/2 

where ^ determines the chemical potential in units of the Fermi energy, /j = ^£/r with £/r = 
kp/{2m). The two NLO parameters ci,C2 are related to the momentum dependence of correla- 



tion functions. The phonon dispersion relation, for example, is given by llioll 



qo = Vsq 



,2 /^/„, , 3^ \ r , ^/ 4, „/ 2^ 



1 -7rV2^ ( CI + ^C2 ) ^ + 0(^^log(^^)) 



(4) 



where v.^ = y ^/3vf is the speed of sound and vf = kf jm is the Fermi velocity. The static suscep- 
tibility is defined by 



-li^q) = -i I dtd^xe^'^'^{\\f'^\\f{0)\\f'^\\f{t,x)) 



(5) 



where \|/' \|/ = V|/aV" (a = 1 , 2) is the sum of the spin up and down densities. The susceptibility is 



related to a different linear combination of ci and C2 lIlOll 



X(^) 



mkf 



1+2712^/2^ ( CI - ^C2^ ^ + 0{q'\og{q^)) 



(6) 



The effective lagrangian can be used to compute the groundstate energy of fermions confined by 
an external potential. The energy of A/^ fermions in a spherically symmetric trap V{x) = jinoy^x^ is 



E = ^^a)(3A^)4/3 - 3V2n^^(i) fci - ^cz] {3Nf/^ + .... 



(7) 



In this work we will derive an energy functional that depends on the local density n{x). This 
functional is the Legendre transform of the pressure. 



£[n(jc)] =^n{x)-P[iJ~V{x)]. 



(8) 



The energy functional is easily derived from the effective lagrangian. Up to NLO in the derivative 
expansion it is sufficient to consider the tree-level effective lagrangian llioll . The only difficulty 
is to invert the relationship between the density and and the chemical potential. This can be done 
order by order in the derivative expansion. We write 



n\^-V{x)] = no\^-V{x)]+dni\^-V{x)]+d^n2\M-V{x)] + . 
/J-V{x) = iJo[n{x)]+diJi[n{x)]+d'^iJ2[n{x)] + ... 
£[n(jc)] = £o[n(jc)]+5£i[n(x)]-f52£2[n(x)] + ..., 



(9) 
(10) 

(11) 



where 5 is used as an expansion parameter. The functions no,ni,... arise from differentiating 
the leading order, next-to-leading, etc. terms in the effective lagrangian with respect to /j. The 
functions fiQ^/ui,... can be found by inverting this relationship order by order. We find 



£o[n(x)] = n{x)V{x) + iuo[n{x)]n{x) - Po\/uo[n{x)]], 



(12) 



with^o[^W] = ('^o) H"W] ^iid £i[n(x)] = —Pi{iJo[n{x)]). This yields 



£(jc)=77(jc)y(jc) + 



3-22/3 _,/, 4 2ci-9c2(Vn(x))2 12 C2 ^2 



5^/^mcJ 



n{xfl^ 



45 mcQ 



n{x) 



5 mcQ 



V^n(jc). (13) 



The first two terms correspond to the local density approximation (LDA) and the terms propor- 
tional to c\ and C2 are the leading correction to the LDA involving derivatives of the density. We 
note that the last term proportional to V^n{x) does not contribute to the total energy of a finite 
system. 

III. EPSILON EXPANSION 

A. Lagrangian and Feynman rules 



At unitarity the determination of c\ and C2 is a non-perturbative problern, and we will perform 
the calculation using an expansion around J = 4 — £ spatial dimensions Un. \VM- The epsilon 

n 

expansion has proven to be useful in calculating the equation of state [|13|] . the critical temperature 



[|l4|] . few -body scattering observables [1150, and the phase structure of spin-polarized systems [|l6|] 
Our starting point is the lagrangian 



L=^' 



2m 



^+^^^03^+ f^^a+^4) + /i.c.) ^"^^ , (14) 

V /Co 



where ^ = (\j/|,\j;|)^ is a two-component Nambu-Gorkov field, a, are Pauli matrices acting in 
the Nambu-Gorkov space, a± = (oi ±?'a2)/2, (|) is a complex boson field, and Q is a coupling 
constant. In dimensional regularization the fermion-fermion scattering length becomes infinite for 
1/Co^O. 

The epsilon expansion is based on the observation that the fermion-fermion scattering ampli- 
tude near d = A dimensions is saturated by the propagator of a boson with mass 2m. The coupling 
of the boson to pairs of fermions is given by 

(b). We write the boson field as 



m \ 2% 
In the superfiuid phase (|) acquires an expectation value (|)o 
^ = ^o + 8^>- The lagrangian is split into a free part 



^o = 'i' 



V2 



V2 



'^+^ l''^o+4;;;i^' 



(16) 



and an interacting part £/ + L^, where 






4m 



(17) 
(18) 



Note that the leading self energy corrections to the boson propagator generated by the interaction 
term £/ cancel against the counterterms in Let- The chemical potential term for the fermions is 
included in £/ rather than in Lq. This is motivated by the fact that near d = 4 the system reduces 
to a non-interacting Bose gas and /j ^ 0. We will count /j as a quantity of 0(e). The Feynman 
rules are quite simple. The fermion and boson propagators are 



D{po,p) 



Po-Ej 



po + e,: 



-<i}o 



Po 



(19) 
(20) 



Po - ep/2 

where Ep = e^^ + ^q and e^ = p^/{2m). The fermion-boson vertices are igo^. Insertions of the 
chemical potential are ifJO^. Both g^ and /u are corrections of order £. 

We shall make use of the following results that have been obtained at NLO in the epsilon 
expansion [|l2n 

2^ 



<^o 



K 



£ 

1 
£ 

£3/2 



[l + (3C-l+log(2))£ + 0(£^ 



l-i(2Y-l-21og(2)) + 0(£2) 



«4o 

In 



d/2 



l + ^£log(£)-^(12C-5 + 51og(2))£ + 0(£2) 



(21) 
(22) 

(23) 



Here, (|)o is the expectation value of the boson field, n is the density, and ^ determines the chemical 
potential in units of the Fermi energy, /j = ^£/7. The quantity C ~ 0. 14424 is a numerical constant 
that appears in the calculation of the two-loop effective potential, and y ~ 0.57722 is the Euler 
constant. 



B. Phonon propagator 



The phonon dispersion relation at LO in the epsilon expansion was obtained by Nishida in lll4n . 
Here, we briefly review his results. We introduce a two-component scalar field 4> = (cp, cp*). The 
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FIG. 1 : Dyson Schwinger equation for the phonon propagator. Dashed hnes denote the free boson propa- 
gator, double dashed Unes denote the full propagator and n is the boson self energy. An^ows show the order 
in which cp and cp* are contracted. 



scalar propagator is 



-1/ 



'D~\p) = 'D-\p)-n{p) 



[2)-nP)]ll ['^-\p)h2 



(24) 



[2)-l(p)]21 ['D-\p)]22 

where T) is the full propagator, T>o is the free propagator, and H is the self energy, see Fig.[TJ The 
free propagator does not have off-diagonal (anomalous) components. The diagonal terms are 

[V(/')]ii = [2?o^(-/^)]22=Po-|. (25) 

The self energy diagram at LO in the epsilon expansion are shown in Fig. [2l We find 

nil = n22 = -2A/+^ + o(£2), 



n 

At leading order fx = £(|)o/2 and 
The dispersion relation is 



I2 = n2i = ^ + 0(e2). 



po + Y+f^ 



(26) 
(27) 



pl-'iCi+2^) 



-f^ 



-Po + ^+P 



(28) 



Po = -^Jep{ep + 4ij)c^y/^pl\ + ^ + ...] 



(29) 



which shows that the spectrum contains a Goldstone mode with a linear dispersion relation, po 
VsP, where v^ = ^y^/{2m). 
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FIG. 2: Leading order contributions to the boson self energy. Full lines denote fermion propagators in 
the Nambu-Gorkov representation. Arrows indicate the order in which ^ and *P^ are contracted. A cross 
denotes a /^-insertion from £/. There is a contribution from the first counterterm in Let which is not shown 
here. 



C. Static susceptibility 

The one-loop contribution to the static susceptibility (see Fig. [3^) is 

d^k f dko 



i{q) 



-i I -^ I ^Tr[G{k + q/2)a,G{k-q/2)a3] 



(30) 



In J 



d/2 



where we have expanded x{q) in powers of momentum, treating q^ as a quantity of order £. We 
observe that the one-loop contribution scales as %(0) ~ e*^. This should be compared to the thermo- 
dynamic result %(0) = —{dn)/{d^) ~ £^^. In order to get an enhancement by two inverse powers 
of £ we need to consider graphs that contain massless particles. The dominant contribution comes 
from phonons, see Fig. [3h,c. The LO phonon term is 



X{q)=g^n3+{q)'Dn{q)n,^{q)+n,+ {q)'Dniq)n,+ {q)+h.c 



(31) 



where !>, ,■ is the phonon propagator and 



n3±(^) = -il-^I^Tr[G{k + q/2)a,G{k-q/2)Oi] 



(271) 

'S_>„_,o,p„e.i(.-losP))V-l(|-)e|(^ 



(32) 



d/2 



£(^0 I 2' ' '' 8 24 ymd^o^ 

Using the leading order phonon propagator derived in the previous section we find in the static 
limit 



'Dn{q)+'Dn{q) 



1 



>/,4fil)+o(,4« 



2/u + £q/2 2^ 1^ 8 \m/j 



(33) 
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FIG. 3: Fig. a) shows the one-loop contribution to the static susceptibihty. The wavy line denotes an external 
source coupled to \|/^\|/. Figs, b) and c) show the leading order phonon contribution. The double dashed line 
is full phonon propagator defined in Fig. [T] 



We can now determine the static susceptibility 

2 



X(^) = - 



£^ 



l-(Y-log(2))£ + 0(£^) 



8 \mij/ 



2n 



d/2 



(34) 



We can compare the momentum independent term to the prediction from the relation %(0) 
— {dn)/(dij). Using equ. (|22|) and (|2TI) we find, at NLO in the £ expansion, 

d/2 



%(0) = -^|i -^(y-iog(2))£ + o(£2) 

£A/ t 2 



2n 



(35) 



which agrees at leading order, but not at NLO. 



D. Higher order corrections 

The LO phonon dispersion relation and susceptibility depend on 0{e) terms in the boson self 
energy. The NLO phonon dispersion relation requires 0{e^) corrections. Since the LO curvature 
term in the dispersion relation is proportional to p^ / (m/j) we can count po , £p as quantities of order 
£. The one-loop self energies, expanded to NLO in £, are given by 

nil = -(po-|){l-^(Y-log(2))£| + ^|l + i(5-3Y-log(8))£} + ... (36) 
Tin = ^{l + ^(l-y-log(2))£}-^ + ... (37) 



a) ^ b) ^ c) 






FIG. 4: Higher order contributions to the diagonal boson self energy ITii. Fig. a) shows a /^-insertion into 
the one-loop self energy. This diagram is combined with the second counterterm from Let (not shown). 
Figs, b) and c) show the "vertex-type" and "self energy-type" two-loop contributions. The coiTcsponding 
contributions to the off-diagonal self energy YI12 aie not shown. 

and the NLO expression for (|)o is given in equ. (|2TI) . The term — (po — £p/2) is canceled by the 
first counterterm in Let- The ^ insertion into the one-loop self energy is 

nil = -2A/|l-^(l+2y-21og(2))£ + 0(£2)| (38) 

ni2 = ^|l+0(£)|, (39) 

where the term — 2/j is canceled by the second counterterm in Let- There are two two-loop self 
energy diagrams, see Fig. ID We will compute these diagrams in App.lAl The result can be written 
as 

nil = Cl(^oe^ ni2 = C2(^oe^ (40) 

where Ci,2 are numerical constants. These two constants are constrained by some general relations. 
First, the phonon is a Goldstone mode and the dispersion relation has to satisfy ()iq{q = 0) = 0. We 
also know that the velocity of sound is related to the equation of state. In J = 4 — £ dimensions 



T" 



y4(>n+'"'^))- (^'> 



These two conditions determine Ci,2- We find 

Ci = -|c~ -0.64908, C2 = -|c^ -0.21636. (42) 

In App.lAl we demonstrate that these results agree with an explicit calculation of the two-loop self 
energies. With these results, the inverse boson propagator at NLO takes on a very simple form. 
We find 

^-'=Z\'';J~' ^7 , Z=l-l(y-log(2))£ (43) 

-g^-A/ -po-f-/jJ 
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FIG. 5 : Higher order contributions to the static susceptibility. Fig. a) shows the /j-insertion into the one-loop 
diagram, Figs, b) and c) show /j-insertions into the phonon contribution. 



The phonon dispersion relation is 

7^0 = V 



^0+l){ 



! + ■ 



( 



1 



+ . 



(44) 



8a* V 4. 

We now consider NLO corrections to the static susceptibility, see Fig. \5\ The first diagram is 
the /j-insertion into the one-loop graph. This graph is 0(£/j), which is an 0{e'^) correction to the 
leading order term in x(0). The second and third diagram contain /j-insertions into TI3+. These 
diagrams are 0(1), an 0{e^) correction to the leading term. Two-loop corrections to n3_|_ are also 
suppressed by at least two powers of £. This means that, in addition to NLO corrections from the 
one-loop n3_|. already given in equ. (|34l) , NLO corrections to the static susceptibility arise solely 
from higher order corrections to the boson propagator. We find 

2 \ r. -^ / A \ rf/2 



x(^) = - 



ep 



1 



1 



(Y-log(2))£ 



1 



1 

8 \mij 



1 



+ 0(^^) 



2% 



We observe that %(0) matches the NLO prediction from the relation %(0) 
equ. (I35l) . 



(45) 



-{dn)/{dij), see 



IV. MATCHING 

We found that the leading terms in a small momentum expansion of the phonon dispersion 
relation and the static susceptibility satisfy the low energy predictions v^ = {dP)/{dp) and %(0) = 



11 



— {dn)/{dij) at NLO in the epsilon expansion. This means that we can use the curvature terms to 
fix the low energy constants ci and C2- In order to be consistent with the low energy theorems we 
have to perform the matching ind — A — £, dimensions. The low energy effective Lagrangian in d 
dimensions is 



£ = com^/2xi+^/2 + cim^/2-ii^ + . '^ 



(V^cp) -d^mV^V 



Xd/2-\ 



(46) 



The powers of m and X follow from the scaling dimension of the fields. The factor dP' in the C2- 
term is a non-trivial consequence of conformal invariance in d dimensions [110|]. In d dimensions 
the relation between cq and h, = n/Zf is 

2 



Co 



(47) 



(27r)^/2r(2 + J/2)^^/2' 

The two NLO parameters ci,C2 can be related to the momentum dependence of the phonon dis 
persion relation and the static susceptibility. In d dimensions we find 



qo = Vsq 



2 1 



4 / d 

d{d + 2)co \ 2 / m/u 



(48) 



and 



xiq) 



d{d + 2)co ^^^a;2^l/2-i 



! + ■ 



8 



C,-d^^-l]c2] — 



2 



(49) 



d{d + 2)co 

We can now match the curvature terms in equ. (|44|) and P3I) to equ. (|48|) and (|49l ). From the 
phonon dispersion relation we get 



d 



Matching the static susceptibility gives 



d{d + 2)cQ / £ 
64 V ~4 



d ,2 f d 

c\ + -C2 = ci-d I- 



1 C2. 



This implies that cj vanishes to NLO in the epsilon expansion cj/ci = 0{e 
consistent with the constraint C2 > [10]. The ratio ci/co is given by 



^=-^ri-^+... 

CO 8 V 3 



(50) 

(51) 
This is (barely) 

(52) 



At NLO we obtain the following density functional for non-relativistic fermions at infinite scatter- 
ing length 

(53) 



£,{x)=n{x)V{x) + 1.364"^^^ + 0.022^^^^^ + OiV\) 



m 



mn{x) 
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It is interesting to compare this result to the density functional for non-interacting fermions [|l7n 

\2 



££^^(x)=n(x)y(x) + 2.871^^^^ +0.014^^^^ + 0.167 ^^^ + 0(v4n) 



m 



mn[x) 



m 



(54) 



which is known at the "extended Thomas-Fermi model" (ETF). The energy of A'^ fermions in a 
spherically symmetric harmonic trap is 



^(0(3A^)4/3 1 + 



+ ... 



(55) 



4 "^■"■' V' ' (3A^)2/3 
where ^ ^ 0.475 (see equ. (|23l) ) and Cs = — (32ci)/(5co^) ^ 1.68 at NLO in the epsilon expansion. 



The result for free fermions is 



E...^ 1.(3^)*/^ (,+5^ 



+ ... 



(56) 



We observe that the coefficient of the A'^'^' ^ term in the ETF functional is larger than the corre- 
sponding coefficient in the unitarity limit. This simply reflects the fact that the interaction between 
the fermions is attractive and ^ < 1 . What is more surprising is the fact that the ETF functional 
corresponds to a significantly smaller value of Co. Numerical results for uri io N = 30 harmon- 
ically trapped fermions can be found in ni8l . 1191 12011 . For small A'^ the corrections to the local 
density approximation are not very well fit by a N~-^'^ contribution, and the authors of [|l8L Il9|] 
did not attempt to extract ^ and c.v independently. Under the assumption that the data can be de- 



scribed by £■ = ^Eetf they find values ^ 



cepted bulk value ^ ~ (0.40 - 0.44) Ul 



0.47 — 0.50) which are larger than the commonly ac- 



22L 



21 



2411 . On the other hand, accepting the bulk value 
^ = (0.40 — 0.44) implies larger values of c^ than the one predicted by the extended Thomas- 
Fermi model. Taking E{N = 20) = (41.3 - 43.2)(0 from UJ, [ij] and ^ = (0.40 - 0.44) gives 
Cs — (0.9 — 2.5), consistent with our result Cs = 1.68. We note that the data of UM u3, \2m\ are 
even better fit by a functional of the form E{N) = (0\/^/4 ■ {3N)'^/^{\ + c/(3A^) ^Z^). A correction 
of the form \/N^'^ cannot be obtained from a local energy density functional of the form given 
in equ. (fT3l) . nor is it compatible with the structure of non-leading terms in A^ generally assumed 



in the literature 1125 . 



2611 ■ It is an interesting challenge to determine whether more complicated 



functional (see Sec. |V]) can yield corrections to the energy of A'^ harmonically trapped fermions 
that scale as l/N^/'^. 
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V. SUMMARY AND OUTLOOK 

We have computed an energy density functional for dilute non-relativistic fermions at unitarity. 
Our approach is based on an effective field theory of the unitary gas which takes into account the 
effects of spontaneous symmetry breaking, the existence of Goldstone modes, and the constraints 
from Galilean and conformal symmetry. We have used an epsilon expansion to compute the co- 
efficients in the effective lagrangian at NLO in the derivative expansion. Our main result is given 
in equ. (l53l) . It is interesting that at NLO in the epsilon expansion only one of the two possible 
two-derivative terms appears. 

There are several interesting lines of investigation that we wish to pursue in the future. The 
first is the problem of constructing energy density functional that depend on more than one type 
of density. One may consider, for example, the spin density (mostly of interest for applications in 



atomic physics), or the superfluid density 1127 . 



p spm aensity (r 
la y,y]. It i 



way of constructing functionals of the Kohn-Sham type Ml 



is also important to find a systematic 



3211 . or extensions of Kohn-Sham 



theory that contain anomalous densities (as in the Hartree-Fock-Bogoliubov approximation), see 



II27I1 . Finally, it is important to study more realistic interactions for neutron matter, in particular 



the effects of a finite scattering length or non-zero effective range 113 311 . and the effects of explicit 
pion degrees of freedom 
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APPENDIX A: TWO-LOOP SELF ENERGY DL^GRAMS 

In the appendix we compute the two-loop self energy diagrams that contribute to the phonon 
propagator at NLO. The two-loop "vertex-type" diagram shown in Fig. Hb) is given by 

-in'^SiP)=-s'l-^-0y^^[Gn{k)G2MG2^^^ (Al) 

+Gi2{k)Gn{q)G2i{q-p)G22{k-p)D{q-k)], 
M'hp)=-s'l^0^^^[Gn{k)Gn{q)G22{q-p)Gn{k-p)D{k-q) 

+G22{k)Gn{q)Gn{q-p)Gi2{k~p)D{q-k)]. 

In the £-expansion, counting po ~ P^ ~ £ we only need the self-energy correction at zero energy 
and momentum. We get 

-/<)(0)=-/<)(0) (A2) 



7- 



4 f dq^ dk" d d 
ig 



kg + £yt ^0 — Eyt 



{2nY{2nYdqodko [{ko-Ek)^ {qo-Eq)^ko-qo-eq-k/2_ 



/-^a + 0{e'). 



qo ^Eq 

ko -^ —Ek 



The constant a can be determined by performing the integrals in J = 4 spatial dimensions. A 
numerical calculation gives a ~ —0.267359. 

The "self energy-type" two-loop Feynman diagram Fig. H];) is given by 

-/nSf=-2//^^^^[GH(^)GnWG22(^)G22(^-;^)Z)(fc-^) (A3) 

+Gn{q)Gn{k)G2i{k)G22{k-p)D{q-k)]=B + F, 

+G22{k)G2i{k)G2i{k-p)Gxi{q)D{q-k)]=H + K. 
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In the limit p — 0, changing variables ko -^ —ko shows that F — K and F = H. We get 



D 



. 4 f dq" 



{2kY{2kY 2Eq dkl 
( b, 

K 

dq'^ dk"^ Eq-eq 32 



^ J UnY 



{iTzYilnY 2Eq dkl 



iko-EkY{ko-Eq-ek-q/2) 



ih-^M 

{ko + EkYi-ko-Eq-ek-g/2) 



ko = -Et: 



(A4) 



k)=Ek 



.2(^oe2 „ 



Numerical evaluation gives b ~ -0.5822930 and / = /z = fc ~ 0.09742858. Collecting all the 
terms 

n5*^(0) +n[f (0) = ^^^{2a + b + f)=Ci(^oe^ - -0.64908(^oe^ (A5) 



n^7(0) +n^7(0) = ^^^^(a + /) =C2(^oe^ ~ -0.21636(^oe^ 



rWm^^rrW^n^_'^^oe^^_ ^,_^_^_^2^. noi^o^A_.2 



which agree very well with the determination in equ. (1421) . 
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